Single mode heat rectifier: Controlling energy flow between electronic conductors 
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We study heat transfer between conductors, mediated by the excitation of a monomodal harmonic 
oscillator. Using a simple model, we show that the onset of rectification in the system is directly 
related to the nonlinearity of the electron gas dispersion relation. When the metals have strictly 
linear dispersion relation a Landauer type expression for the thermal current holds, symmetric 
with respect to the temperature difference. Rectification becomes prominent when deviations from 
linear dispersion exist, and the fermionic model cannot be mapped into a harmonic- bosonized- 
representation. The effects described here are relevant for understanding radiative heat transfer and 
vibrational energy flow in electrically insulating molecular junctions. 
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Developing nanoscale devices with asymmetric conduc- 
tion properties is a long standing challenge both for fun- 
damental science, for resolving the microscopic mech- 
anisms controlling transport, and for practical appli- 
cations. In the last years, rectification in mesoscopic 
and molecular level systems have attracted much at- 
tention. Recent works demonstrated electrical rectifi- 
cation in various systems, including mesoscopic semi- 
conductor structures one-dimensional (ID) systems 
of interacting electrons Q, and within molecular junc- 
tions [1, Q . The growing interest in the thermal prop- 
erties of nanoscale structures Q turned rectification of 
phononic current into a topic of great interest, with im- 
plications on nanoscale machinery Q and thermal com- 
putation Ml. The different setups demonstrating this ef- 
fect @, I9l.ll0l [TH [l2j all rely on two basic assumptions: 
The device should be asymmetric with respect to the two 
terminals, and the system's normal modes should nonlin- 
ear ly interact. 

Radiative thermal conductance, where heat exchange 
between metals is mediated by the generation of photons 
is a new topic of interest [l3[ . It was recently proved that 
at low temperatures photonic heat conduction is quan- 
tized [l4j |. setting an upper bound on single-channel in- 
formation flow [15| . From the practical aspect, at the 
nanoscale, radiative heat flow may compete with vibra- 
tional energy transfer, thus these two processes must be 
considered for properly estimating the thermal conduc- 
tance of molecular level systems [13| . 

In this paper we investigate monomode mediated en- 
ergy exchange between two metals, and resolve the nec- 
essary conditions for manifesting thermal rectification. 
We show that when the metallic leads have linear disper- 
sion relation, or in other words, when the bosonization 
approach can be employed to yield a harmonic Hamil- 
tonian, a Landauer type expression for the energy cur- 
rent holds, symmetric with respect to the temperature 



difference. In contrast, for conductors with nonlinear 
dispersion relation, when the approximations involving 
the standard bosonization scheme break, the energy cur- 
rent can be rectified. Deviations from the Tomonaga- 
Luttinger bosons picture [l6| thus eminently relates to 
the onset of rectification in the system. 
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FIG. 1: A schematic representation our model. While electron 
flow is blocked, even for > Hr, heat current is flowing right 
to left (Tr > Tl) through an excitation of the local harmonic 
mode. The curved arrows represent energy transfer processes 
between the leads and the intermediate mode. 

The fermionic model consists two metallic leads held 
at different temperatures, coupled by a single harmonic 
mode. The Hamiltonian includes three contributions 
Hp = Hq + H e + Vf, 
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Hq includes a single harmonic mode (subsystem) of fre- 
quency u>o, creation operator &J. It can be considered as 
a monomodal electromagnetic field, or in a different con- 
text, it represents a local vibrational mode of an interme- 
diate molecular unit. Henceforth we refer to it as a local 
mode. This mode interacts with two fermionic reservoirs 



2 



(H e ), where cj, k (c Vtk ) creates (annihilates) an electron at 
the v = L, R metal with momentum fc, disregarding the 
spin degree of freedom. The oscillator-metals interaction 
term Vp couples scattering processes within each metal 
to the subsystem degrees of freedom, where / is a local 
mode operator. For simplicity, we use / = bg+bo, and as- 
sume that the coupling constants a v are energy indepen- 
dent and real numbers. Since Vp does not directly couple 
the reservoirs, the model ([1]) leaves out charge transfer 
processes, assuming the tunneling barrier is high and the 
constriction is long. However, energy can be transferred 
between the two metals mediated by the excitation of the 
local mode. For a schematic representation see Fig. [TJ 

While the aim of this paper is fundamental, to under- 
stand the microscopic origin of asymmetric transport in 
simple fermionic Hamiltonians, the generic model stud- 
ied here may be realized in different systems: (i) Ra- 
diative energy transfer, where the transfer of the lead's 
excess energy is mediated by the excitation of electromag- 
netic modes. Other mechanisms for energy transfer are 
neglected, e.g., thermionic emission, vibrational energy 
flow, and photon tunneling, (ii) Vibrational energy flow, 
where thermal energy is transferred between two metal- 
lic terminals through a vibrating link. Here the basic 
assumptions are that the molecular link electrically in- 
sulates, and that vibrational energy cannot directly flow 
from the leads to the molecule due to a large mismatch 
of the metal-molecule phononic spectra. 

For ID noninteracting electrons with an unbounded 
strictly linear dispersion relation, e k = ep + vp (|fc| — kp ) , 
tp is the Fermi energy and vp is the velocity at the Fermi 
energy, bosonization can be employed to yield an equiva- 
lent bosonic Hamiltonian, H B — Hq +-H& + Vb- Here Hq 
is the same as above, while the fermionic degrees of free- 
dom are written in terms of Tomonaga-Luttinger (TL) 
bosons [l6j, creation operator b\ (y — L,R), 

V B =J2K„,M, g + K, q )(bi + b ). (2) 



The new coup ling parameters relate to the couplings 
a v [Eq. (TT])] [13, Assuming an ohmic dissipation, 
the spectral function of the v boson bath g v (uS) can be 
written in terms of the fermionic parameters, 

g„{u>) = 47ry\y i<7 fl(c<; - cj q ) = 2ttuj^; 
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with p u {eF) the density of states at the Fermi energy. 
The boson Hamiltonian Hb is fully harmonic, and thus 
can be written in terms of noninteracting collective 
modes. The heat current in this model can be exactly 



calculated to yield a Landauer type expression 

J = - J T{uj)[n L B {uj) - n%(w)]wdu). (4) 

Here n B (u>) = [e"/ T " — 1] _1 is the Bose-Einstein occupa- 
tion factor with temperature T v . The L to R transmis- 
sion coefficient is given by T(w) 
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[(u>a-wg)*+(r£+rg)^] 
[19J, with the relaxation rate T B = 2tt J2 q K t q^( w ~ w g)- 
In the weak coupling limit, T B < uiq, the transmission 
coefficient is sharply peaked around ujq, and Eq. Q re- 
duces into a resonant energy transfer expression, 
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Here T B is calculated at the (local oscillator) frequency 
uiq. For weak coupling, irp(ep)a < 1, Eq. yields 
£ = 2p(e F ) 2 a 2 , and the relaxation rate is given explicitly 
by 

T B = 27rcu alpUe F ). (6) 

Eqs. (JU-JU describe the heat current of the model 
Hamiltonian ([1]) under the assumption that the two 
fermionic reservoirs have a strictly linear dispersion re- 
lation. It is clear that in this case thermal rectification 
cannot show up, since exchanging the temperatures of 
the leads simply reverses the sign of the heat current Q , 
not the absolute value [2(| ■ We show next that when the 
bosonization method cannot be trivially employed, i.e. 
when the dispersion relation is not linear, the model ([T]) 
can bring in an interesting rectifying behavior. 

Assuming a general dispersion relation, the dynamics 
of the model |T]) is analyzed in the weak coupling system- 
bath limit. Going into the Markovian limit, the proba- 
bilities P n to occupy the \n) state of the local oscillator 
satisfy the master equation 



(7) 



where the transition rate from the local oscillator state 
\m) to \n) is additive in the L and R reservoirs, k n —>m = 
^n— >m + kfL^fn ' due to the linear form of the interaction 
[Eq. (JTJ>] 10]. In steady state, the heat current across 
the system is given by (calculated e.g. at the L side), 
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formula, the transition rates are given by 

K^ m = 27r|/„ Mi | 2 ^ \a v> k; V ,k>\ 2 
fc.fc' 

xn F (e k )[l - n F (c k >)]5(e k - e k > - E m . n ) 

= 2^|/„ Mi | 2 J den F (e)[l - n v F (e - E m , n )]F u (e). 



-2n\f m ^ n \ 2 n B (E m . n ) 
x / de[n F (e)-n F (e-E m>n )]F v (e). 



(9) 
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Focusing on the last equality, interestingly we see that the 
thermal properties of the reservoirs are concealed both 
within the Fermi-Dirac distribution function n F (e) = 
[ e O-Ai.>)/T„ -|_ i]- 1 and the Bose-Einstein occupation fac- 
tor n B (e) = [e e / T,y — It is therefore clear that when 
the integral yields a temperature independent constant, 
the statistic of the reservoirs is fully bosonic. The other 
elements in ([9]) are the matrix elements of the system op- 
erator f m<n = (m\f\n), and the dimensionless interaction 
term F„(e) = \a v \ 2 p^(e)p u (e — E m n ), which encloses the 
properties of the reservoirs multiplied by the system-bath 
couplings a v . These couplings might be taken different 
at the two ends, bringing in a spatial asymmetry. 

We assume next that the density of states slowly varies 
in the energy window E m ^ n . The interaction function is 
then expanded around the chemical potential (2lj . 
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with 7„ a dimensionless number of order unity. Using this 
form, the integration in Eq. ([9]) can be performed when 
the Fermi energies are much bigger than the conduction 
band edge, p v 3> E c [2l|. Making use of the following 
relationships, 



de[n F (e) - n F (e - E m>n )] « 
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we get 



F v {p v ) + lAj u ^ 

Pv_ 



(12) 



Note that n B (—E„ ltn ) = ~[n B (E mn ) + 1], thus the ex- 
citation and relaxation rates induced by the v reser- 
voir satisfy the detailed balance relation, k^ m /k^ n ^ tn — 
g-Em.n/Tv _ consider next the limit of a constant den- 
sity of states, taking 7 = 0. Eq. (fT"2|) then becomes 

K->n-i = nr^(w )[l + n B (uJo)}; 

T f (lu ) = 2-nF v {p v )u , (13) 

where the bilinear interaction form was employed / = 
(6q + 60), leading to nearest neighbors transitions only. 
We can also calculate the states population in steady 
state by putting P n = in Eq. 



x n (l-x); x 
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Finally, the heat current for the 7 = case is calculated 
with the help of Eq. (JSJ, 



J = LU 



-iJ^[n|M-n*M]. (15) 
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T F is calculated at the frequency ujq. In the linear disper- 
sion limit we thus recover the resonant energy behavior 
([5|) obtained with the equivalent boson Hamiltonian 
Note that the rates calculated with the different methods, 
T F and T B , are equal, see Eqs. ([6]) and (fT3|) . 

We evaluate next the current when 7 7^ 0, i.e. for a 
model with an energy dependent density of states. For 



the linear coupling case, / 



6q, only transitions 



between nearest levels survive and Eq. 



fT2]) becomes 

K-+n-i = nT F {w )[l + n B {u )] (l + \ u fy, where \ v = 
1 ' 1 ? The steady state population is therefore given by 



Kin*)' 

Eq. dHJ) with Y F -> T F (l+\ v T v /p v ). Next we calculate 
the energy current using Eq. ©. In the classical limit 
(T v > ujq) we obtain 



J = 



T F T§(1 + X 



A B ^' 



r^(i + A i ^) + rf(i + A fl ,^) 



(T L -T R ). (16) 



The dimensionless parameter X v effectively quantifies the 
deviation from the linear dispersion case. System-bath 
interactions a v are contained within the coefficients T F 
computed at the subsystem frequency luq. 

Comparing Eq. (TIB)) to the high temperature limit of 
Eq. (|15p interestingly reflects the deviations from the 
TL boson picture in transport properties. Eq. (|16|) has 
three important characteristics: First, the heat current 
obtained is a nonlinear function of the temperature differ- 
ence. In fact, for T F > Pp, J oc aiAT + Xa 2 AT 2 . The 
constants a,\^ depend on subsystem and bath parame- 
ters, AT = Tl — Tr. Second, this expression manifests 
nontrivial controlability over the energy current by tun- 
ing the chemical potentials. Thirdly, this result demon- 
strates thermal rectification, as the thermal current is 
different when switching the temperature bias, assuming 
some asymmetry is included, e.g. the chemical poten- 
tials, or the local mode-bath interactions are different 
at the two terminals. When the electronic properties 
of the reservoirs are equivalent, A = A„, p a = Pv, we 
find A J = J+at + J- at oc -A^(p£ - r£), where 
J±at = J(T L - T R = ±AT). The proportionality factor 
is given by G w T F T^/{r L F +T^) 2 for XT v /p a < 1. Thus, 
a system made of a local harmonic mode coupled asym- 
metrically to two metals with energy dependent density 
of states, rectifies heat. Rectification becomes more ef- 
fective with increasing A, i.e. when the density of states 
strongly varies with energy. Therefore, the same param- 
eter that " quantifies" the departure of the fermion model 
from the boson picrture, determines the effectiveness of 
rectification in the system. Figure [5] presents this effect. 
We plot the absolute value of the ratio J+at/ J- at, and 
find that it linearly departs from unity with increasing 
temperature difference and A, in agreement with the an- 
alytical expression (T a = Tl -^ Ta ), 
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FIG. 2: Thermal rectification with increasing deviations from 
constant density of states. A = (full); A = 0.1 (dashed); A = 
0.5 (dashed-dotted); A = 1 (dotted). Y% = 0.5, T§ = 0.02, 
l^a — /ii/ = 1, temperature of the cold bath is Ti ow =0.1. 
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FIG. 3: Thermal rectification due to chemical potential dif- 
ference. AT = 0.23 (full); AT = 0.15 (dashed); AT = 0.01 
(dotted). r£ = r£ = 0.1, fj, L = 1, Ai = Ar=1. 



Figure [3] demonstrates the dependence of the rectifying 
behavior on the leads' chemical potentials for a system 
symmetrically coupled to the terminals using Eq. (fT())) . 
We find that at small temperature difference rectification 
is negligible, while for AT//i£ > 0.1 the asymmetry in 
the current can be large up to 10%. We note however 
that tunneling of electrons becomes a significant source 
of noise at large voltage bias and/or temperature dif- 
ferences. Interestingly, the energy current can be also 
controlled by tuning A/i, an outcome of the energy de- 
pendence of the density of states. 



The thermal current can be calculated for other 
models besides (fTUl) . expressing Eq. © as fc^_> OT = 
-^Afm,n\ 2 n v B {E m ^ n )alg v (T v ), where g„{T v ) is defined 
through this relation. Eq. (fTT|) then reduces to 1Z = 
\J+at/ J~at\ « g R (T R )/g R (T L ) for the linear coupling 
model if ql ^> a R . We consider a superconducting R 
lead (L can be a normal metal), as it effectively sup - 
presses normal electronic thermal conductance [5J, Il4j . 
and numerically evaluate g R (T) with the density of states 

<°«( £ ) = 7P=s5 Q (l e l - A )- Usin S A = °- 2 meV ( alu " 

minum), ujq = 0.4 meV, and T = 0.08 — 0.2 meV, we get 
gn(T) ~ a\ + CL2T where 3> a\. The rectification 
ratio then becomes 7Z ~ Tr/Tl, which can be as large 
as 2.5 for the above parameters, potentially measurable 
using the setup of Ref. [22| . Similarly, a narrow bandgap 
(0.1 eV) semiconducting lead yields 72. ~ 10 for uj = 0.1 
eV at the temperature range T ~ 0.01 — 0.5 eV. 

Previous studies on phononic heat transfer in ID 
chains demonstrated that anharmonic interactions are 
crucial for manifesting thermal rectification [f| 0, [l(| [H| • 
Our result complies with this observation: A fermionic 
reservoir with a nonlinear dispersion relation can be rep- 
resented by a bath of bosons comprising of nonlinear in- 
teractions [23]. The fermionic Hamiltonian |T]) can thus 
be mapped into a fully bosonic model describing a har- 
monic link bilinearly connected to anharmonic thermal 
baths. 

In summary, we present here a simple model of single- 
mode energy transfer between metals. The model can 
describe energy flow through vibrating link and radia- 
tive heat transfer. We resolve the microscopic require- 
ments for manifesting rectification: For noninteracting 
electrons, rectification appears when the reservoirs' den- 
sity of states, or analogously, system-bath couplings, are 
energy dependent. The effects discussed here are not 
limited to the specific form used (fTU)) . and are a general 
manifestation of the breakdown of the TL boson picture. 
Since dissipative reservoirs typically contain anharmonic 
interactions, finite rectification of the energy current be- 
tween metals and dielectric surfaces is an inevitable ef- 
fect. 
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